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ABSTRACT 

We present an analysis of fine-structure transitions of Fe II and Si II detected in a high-resolution 
optical spectrum of the afterglow of GRB 051111 (z = 1.54948). The fine-structure absorption features 
arising from Fe II* to Fe II****, as well as Si II*, are confined to a narrow velocity structure extending 
over ±30 km s _1 , which we interpret as the burst local environment, most likely a star forming region. 
We investigate two scenarios for the excitation of the fine-structure levels by collisions with electrons 
and radiative pumping by an infra-red or ultra-violet radiation field produced by intense star formation 
in the GRB environment, or by the GRB afterglow itself. We find that the conditions required for 
collisional excitation of Fe II fine-structure states cannot be easily reconciled with the relatively weak 
Si II* absorption. Radiative pumping by either IR or UV emission requires > 10 3 massive hot OB 
stars within a compact star-forming region a few pc in size, and in the case of IR pumping a large dust 
content. On the other hand, it is possible that the GRB itself provides the source of IR and/or UV 
radiation, in which case we estimate that the excitation takes place at a distance of ~ 10 — 20 pc from 
the burst. Detailed radiative transfer calculations are required in order to verify that excitation of the 
low-ionization fine-structure states is possible given the intense UV flux from the burst. Still, it is clear 
that GRB absorption spectroscopy can provide direct information on the mode and conditions of star 
formation at high redshift. 

Subject headings: gamma-rays:bursts — ISM:abundances — ISM:kinematics 


1. INTRODUCTION 

Understanding the physical conditions in star forming 
regions at high redshift plays a crucial role in our ability 
to trace the evolution of star formation and the associated 
production of metals. Unfortunately, it is exceedingly dif¬ 
ficult to probe these regions using background quasars be¬ 
cause of their compact size, and possible dust extinction. 
Absorption studies of bright star forming galaxies (e.g., 
the lensed Lyman break galaxy MS 1512—cB58; Pettini 
et al. 2002) provide greater promise because they probe 
the regions where star formation is taking place, but even 
in these rare cases information is only available on the inte¬ 
grated properties and individual regions cannot be probed. 

Gamma-ray bursts (GRBs), on the other hand, are 
known to be the end product of massive stars, and as a 
result their bright afterglows can be used to trace the inter¬ 
stellar medium of their host galaxies. More importantly, 
since GRB progenitors are short-lived, they are likely to be 
buried within star-forming regions undergoing active star 
formation (Bloom et al. 2002). As a result, GRBs may pro¬ 
vide an ideal direct probe of the physical conditions in high 
redshift star forming regions. Over the last several years, 
and particularly with the advent of the Swift satellite, sev¬ 


eral GRB absorption spectra have been obtained revealing 
a large fraction of objects with unusually large neutral hy¬ 
drogen column densities even for damped Lya (DLA) sys¬ 
tems, N 10 22 cm -2 (e.g., Vreeswijk et al. 2004; Berger 
et al. 2005). In addition, it has been noted in several cases 
that the column densities of non-refractory elements (e.g., 
Zn) are higher than in QSO-DLAs, with a depletion pat¬ 
tern suggestive of a large dust content (Savaglio & Fall 
2004). These properties suggest that at least some GRB- 
DLAs are physically different from QSO-DLAs, and that 
they may represent individual star-forming regions with at 
least a modest dust obscuration. 

Given the possible association with star-forming regions, 
several authors have also investigated the impact of the 
GRB prompt and afterglow emission on the local envi¬ 
ronment, particularly in the context of dust destruction, 
photo-ionization, and dissociation and excitation of #2 
molecules (Draine & Hao 2002). The basic conclusion of 
these studies is that bright afterglows are capable of sig¬ 
nificantly modifying their local environment on a scale of 
about 10 19 cm. The detection of such effects can therefore 
provide a unique diagnostic of the immediate environment 
of the burst, and in addition may address the issue of rel- 
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atively low dust extinction in GRB optical afterglows. 

In this Letter we present a high-resolution absorption 
spectrum of GRB 051111, which reveals strong absorp¬ 
tion features from Fe II and Si II fine-structure states 
at a redshift 2 = 1.54948. This is the first example of 
an extra-galactic sight line that reveals the full range of 
Fe II ground-level fine-structure transitions, suggesting 
that their excitation is intimately related to the burst en¬ 
vironment or the burst itself. We investigate both possi¬ 
bilities and their implications in the following sections. 

2. OBSERVATIONS 

GRB 051111 was detected by Swift on 2005 November 
11.250 UT (H. Krimm et al. 2005), and the optical af¬ 
terglow was discovered with the ROTSE-IIIb robotic tele¬ 
scope 27 s after the burst with an unfiltered magnitude of 
about 13 (Rujopakarn et al. 2005). Follow-up observations 
in the first hour after the burst indicate an optical flux 
evolution, F„(R) « 280(f/l hr) -0 9 fiJy (Smith & Swan 
2005). 

Spectroscopic observations of GRB 051111 were initi¬ 
ated approximately 1 hr after the burst, using the High 
Resolution Echelle Spectrometer (HIRES) mounted on the 
Keck I 10-rn telescope (Hill et al. 2005). A total of 5400 
s were obtained in three exposures using a 0.86" wide 
slit. The wavelength range is 4200 — 8400 A. The spectra 
were reduced using the Makee pipeline routines (Version 
4.0.1 of May 2005), which includes optimal extraction of 
orders, sky subtraction, and wavelength calibration from 
Th-Ar arc lamp exposures, including a heliocentric veloc¬ 
ity correction. The orders within individual frames were 
traced using a median combined total of the exposures, 
and atmospheric absorption features were removed with 
the Makee pipeline. A final resampling of the spectrum 
and continuum fitting was performed using the IRAF task 
continuum. 

The spectrum reveals strong absorption features at a 
redshift, z = 1.54948 ± 0.00001, which we interpret to 
arise in the host galaxy (see also Prochaska 2005 8 ). A full 
analysis of the host galaxy and intervening absorption sys¬ 
tems is presented in a companion paper (Penprase et al. 
2005). Here we focus on the detection of Fe II and Si II 
fine-structure transitions. 

3. FE II AND SI II FINE-STRUCTURE TRANSITIONS 

In Figures 1-5 we plot the absorption profiles of all 
transitions of the Fe II ground and fine-structure states 
as a function of velocity relative to the systemic redshift 
of z = 1.54948. The line identifications, observed wave¬ 
lengths, and equivalent widths are summarized in Table 1. 
Also included in the Table are the absorption lines from 
the ground-state and fine-strcuture level of Si II. We de¬ 
tect absorption features from all four fine-structure levels 
of the Fe II 3 d 6 4s 6 D level, (Fe II* to Fe II****), corre- 
spoding to J-values of 7/2, 5/2, 3/2, and 1/2 (see also 
Prochaska 2005). The absorption profiles exhibit a simple 
and symmetric velocity structure extending from about 


—30 to +30 km s -1 and centered on the systemic veloc¬ 
ity of the burst. Strong lines of the Fe II ground-state 
level exhibit an overall velocity range of about —150 to 
+50 km s -1 , but about 85% of the gas is contained in the 
narrow velocity component (Prochaska 2005, Penprase et 
al. 2005). We interpret the dominant component to arise 
from the local environment of the burst based both on the 
kinematic coincidence and the unusual detection of fine- 
strcuture transitions. 

The detection of the full range of Fe II ground-level 
fine-structure states is highly unusual, and to our knowl¬ 
edge represents the first such example in an extra-galactic 
source 9 . In the Milky Way, such transitions have been 
detected in the context of a few dense environments such 
as in the circumstellar disks around (3 Pictoris (Kondo & 
Bruhweiler 1985) and 2 Andromedae (Cheng et al. 1997), 
and the HH 47A bow shock (Hartigan et al. 1999). Ab¬ 
sorption from Si II*, on the other hand, is more prevalent 
and has been previously detected in several GRB-DLAs 
(Savaglio & Fall 2004; Vreeswijk et al. 2004; Berger et al. 
2005), as well as in the lensed Lyman break galaxy MS 
1512—cB58 (Pettini et al. 2002), but not in any QSO- 
DLAs (Howk et al. 2005). This suggests that Si II* may 
trace regions of intense star formation, which tend to be 
missed in quasar sight lines due to the cross-section selec¬ 
tion effect or possibly the associated dust extinction. 

We measure the column densities of the various levels 
using the apparent optical depth method (Savage & Sem- 
baclr 1991). This method has the advantage that it makes 
no a priori assumptions about the functional form of the 
velocity distribution, and at the same time it incorporates 
information from a wide variety of lines with different oscil¬ 
lator strengths, allowing a robust measure of the column 
density in the line cores. The column densities derived 
from this method by integrating over the velocity range 
of ±30 km s _1 , as well as the ratios of the fine-structure 
levels compared to the ground state, are summarized in 
Table 2. 

3.1. Collisional Excitation 

The excitation of the fine-structure levels can be 
achieved either by collisions with ambient electrons, or by 
radiative pumping due to a local radiation field. Colli¬ 
sional excitation of Fe II into its fine-structure states has 
been studied by Keenan et al. (1988), and more recently 
by Silva & Viegas (2002); the latter authors also inves¬ 
tigate the conditions for Si II* collisional excitation. In 
Figure 6 we plot contours of the observed ratios of Fe II 
fine-structure columns relative to the ground-state in the 
T — n e space, and find that the required electron density is 
at least 2 x 10 3 cm -3 . Similarly, the detection of Fe II**** 
lines requires a temperature of at least 400 K. Taking into 
account all four ratios we find that the best-fit density and 
temperature are logn e « 3.5 — 4 and logT « 2.7 —3.1. We 
note, however, that the ratio of Fe II* to Fe II is about a 
factor of two lower than expected based on the ratios of 
the three other excited states. 

A more significant discrepancy is observed in the ratio 


8 A GRB Coordinates Network Circular by Prochaska (2005) provides a cursory list of some properties of the host galaxy absorption system, 
including the detection of fine-structure lines. 
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of the Si II* column relative to that of the ground state, 
AT(Sill*)/iV(Sill) jC 0.06. This is an upper limit since the 
single Si IIA1808 line detected in the spectrum is saturated. 
A comparison to other absorption features from which we 
can determine the Doppler parameter b (Penprase et al. 
2005) suggests that the ratio is about 0.022. Using the 
analysis of Silva & Viegas (2002) we find that for the tem- 
pearure range of 10 3 — 10 4 K, the observed ratio indicates 
an electron density of only <; 20 cm -3 . Conversely, using 
the electron density inferred from the Fe II transitions we 
find an expected column density ratio for Si II* of about 
0.5 — 1, at least an order of magnitude larger than detected. 

The significant difference in physical conditions derived 
from the Fe II and Si II fine-strcuture absorption suggests 
that they either arise in different regions of the absorb¬ 
ing cloud, or that a different mechanism, namely radiative 
pumping is at play. We consider the former possibility 
unlikely since there is no clear reason why Fe II and Si II 
should be segregated within the absorber, particularly in 
light of their similar velocity structure. 


3.2. Excitation by Radiative Pumping 

In the context of radiative pumping two possibilities ex¬ 
ist. First, the fine-structure states can be directly ex¬ 
cited by infra-red radiation. Successive levels starting 
with Fe II* require energies of 384.79, 667.68, 862.613, 
and 977.053 cm -1 , respectively, for excitation. The Si II* 
level requires an energy of 287.24 cm -1 . In the case of IR 
pumping the ratio of the fine-structure states relative to 
the ground state is given by (e.g., Srianand & Petitjean 
2001 ), 

N(X*) 2 n x 

N(X) 1 + n A ’ U 


where the number of photons at the appropriate excitation 
energy is, 


_ I v A 3 
UX ~ 8nhc' 


( 2 ) 


Based on the observed column density ratios listed in 
Table 2 and the energies of the various excitation levels 
we infer the following specific intensities at the location of 
the Si II and Fe II absorbers: 


148.62 x 10 12 Hz) £ 3.9 x 10" 9 erg cm" 2 s” 1 Hz" 1 (3) 

141.15 x 10 13 Hz) « 7.0 x 10" 9 erg cm -2 s _1 Hz -1 (4) 

I v { 2.00 x 10 13 Hz) « 2.2 x 10" 8 erg cm" 2 s' 1 Hz" 1 (5) 

/42.59 x 10 13 Hz) « 3.3 x 10" 8 erg cm -2 s _1 Hz -1 (6) 

142.93 x 10 13 Hz) « 2.2 x 10" 8 erg cm- 2 s -1 Hz" 1 (7) 

These values roughly define a spectrum F v oc v 22 with a 
possible turn-over at v$ ;> 6.4 x 10 13 Hz. 

A second possibility is indirect UV pumping in 
which decays from short-lived excited UV levels popu¬ 
late the fine-structure levels in addition to the ground 
level. Flannery et al. (1980) investigate this process 
for the Si II ion and find that UV pumping domi¬ 
nates over collisional excitation as long as 10 rj = 8.5 x 


10- 14 ?z e (T , e /10 4 K)- 1 / 2 (i?p C /ri7) 2 /- 1 < 1. A detailed in¬ 
vestigation of this process for Fe II has not been performed 
to date. 

4. STAR FORMATION VERSUS GRB EXCITATION OF THE 
FINE-STRUCTURE LEVELS 

As discussed above, the mis-match in electron densities 
required to explain the observerd ratios of Si II and Fe II 
fine-strcuture levels relative to the ground state suggest 
that collisional excitation is unlikely to be the source of 
the observed fine-structure population. In addition, we 
expect that if the excitation is due to electron collisions 
produced as a result of the GRB shock wave, then the 
velocity widths of the fine-structure lines should be sig¬ 
nificantly broader than the observed ±30 km s _1 . Thus, 
excitation due to radiative pumping appears to be a more 
likely scenario. In this context it is instructive to consider 
two sources for the IR or UV photons: star formation and 
the GRB afterglow itself. 

In the context of IR emission from stars, the inferred 
spectrum F v oc v~ 2 ' 2 with a possible turn-over at v$ ;> 
6.4 x 10 13 Hz may be indicative of a thermal spectrum with 
a peak of about 640 K (or greater if the turnover is not 
real), arising from dust re-radiation of absorbed star light. 
The much higher temperature compared to, for example, 
the IR emission from Arp 220 with T « 30 K, as well as the 
required large luminosity, L„s:4x 10 3O (r/lpc) 2 erg s _1 
Hz- 1 (approaching 10 10 Lq), indicate an intense radiation 
field in the environment of GRB 051111 compared to those 
in local star-forming regions. The star formation rate re¬ 
quired to produce the inferred luminosity is in excess of 
100 M 0 yr _1 (Kennicutt 1998), and the high temperature 
requires a large fraction of hot OB stars and a large dust 
content. Perhaps a similar analogue is the super star clus¬ 
ter in the dwarf galaxy NGC5253 (L ~ 10 9 L 0 ), which 
contains ~ 5000 massive O stars within a region of ~ 1 
pc. We note that if this is the correct scenario, then it is 
likely that the GRB had to destroy dust along the line of 
sight in order for the optical afterglow to be detected on 
Earth. 

In the context of UV radiative pumping the condition 
77 < 1 for Si II* excitation described in §3.2 may be satis¬ 
fied in the vicinity (~ 0.1 pc) of OB stars (Flannery et al. 
1980). Thus, a large concentration of such stars may lead 
to the observed Si II* column density. The conditions for 
UV pumping of Fe II fine-structure states have not been 
investigated previously and are beyond the scope of this 
paper, but we expect that a more intense UV radiation 
field is required for excitation of Fe II**** compared to 
Si II*. Given that UV and IR pumping require similar 
star formation conditions, it is likely that both mecha¬ 
nisms operate within the star-forming region. 

Finally, it is possible that the radiative pumping is due 
to the GRB afterglow emission. In this context we use the 
observed flux of the afterglow, F„ ~ 1 —10 mJy in the first 
hour, along with the IR fluxes inferred in §3.2, to find that 
the burst-absorber distance is about 10 — 20 pc, which is 
reasonable for the size of a star-forming region. On the 


1 "The expression assumes a flat spectrum which may not be directly applicable if the UV pumping is due to the GRB UV flash or afterglow 
emission, but we expect the difference to be minor. 
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other hand, it is difficult to reconcile the observed shape 
of the spectrum (F v oc v 2 ) with that of a GRB afterglow 
{F v oc v 1 / 3 ), unless the excitation takes place in different 
layers within the absorber. In the case of UV pumping 
by the GRB afterglow we find that for the observed flux 
this process dominates over collisional excitation within a 
region of about 10 pc unless the electron density exceeds 
10 9 cm“ 3 . 

One possibility in the case of excitation by the GRB and 
afterglow radiation is that the line strengths will change 
as a function of time. We repeated the column density 
calculation for each individual exposure and found no sig¬ 
nificant evolution from t « 1.2 to 2.2 hours after the burst. 

5. CONCLUSIONS 

The high-resolution spectrum of GRB 051111 reveals 
a wide range of absorption features from fine-structure 
states of Fe II and Si II at z = 1.54948. We associate the 
absorption with the local environment of the burst based 
on the fact that the excitation of these fine-structure states 
requires large densities or intense IR/UV radiation fields, 
which are not common in typical interstellar environments. 
This conclusion is also supported by the kinematic coinci¬ 
dence with the burst systemic redsliift. 

A comparison of the conditions required to excite the 
Fe II and Si II fine-structure levels indicates that collisional 
excitation is not likely to be the dominant mechanism, un¬ 
less the two ions are segregated in different regions of the 
absorber. A more likely mechanuism is direct IR pump¬ 
ing or indirect UV pumping. In this scenario we find that 


if the radiation field is due to star formation activity in 
the burst environment, then this requires a large concen¬ 
tration of several thousand OB stars in a compact region 
(~ few pc), perhaps reminiscent of the super star cluster 
in the dwarf galaxy NGC5253. If the UV and/or IR ra¬ 
diation are supplied by the GRB, on the other hand, then 
the distance to the absorber is likely to be about 10 — 20 
pc. This is consistent with calculations by Draine & Hao 
(2002), which indicate that within ~ 3 pc, the burst will 
ionize its environment, leading to a very low column of 
low-ionization states such as Fe II and Si II; outside of this 
region it is likely that Fe II and Si II survive and may be 
excited to fine-structure levels instead. 

Clearly, detailed radiative transfer calculations, along 
with refined calculations of Fe II fine-structure excitation, 
are required in order to distinguish between the different 
scenarios. Still, regardless of the exact details it is clear 
that GRB afterglows allow us to directly probe the con¬ 
ditions within individual star forming regions, and may 
therefore provide direct information on the mode of star 
formation across a wide redshift range. The absorption 
spectrum of GRB 051111 seems to indicate that compact, 
dusty, and dense star forming regions may be prevalent, 
at least in the context of GRB progenitor formation. 

We thank Edward Jenkins and Bruce Draine for help¬ 
ful discussion. EB is supported is supported by NASA 
through Hubble Fellowship grant HST-01171.01 awarded 
by the Space Telescope Science Institute, which is operated 
by AURA, Inc., for NASA under contract NAS 5-26555. 
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Table 1 

Line Identifications for Fe II and Si II 


•^obs 

(A) 

Line 

(A) 

fij 

W 0 

(A) 

4609.49 

Si II 1808.0130 

0.00219 

0.22796 

4632.22 

Sill* 1816.9285 

0.00166 

0.03723 

5736.02 

Fe II 2249.8768 

0.00182 

0.09490 

5763.81 

Fe II 2260.7805 

0.00244 

0.13672 

5935.47 

Fe II** 2328.1112 

0.03450 

0.07175 

5949.25 

Fe II* 2333.5156 

0.07780 

0.15565 

5962.53 

Fe II*** 2338.7248 

0.08970 

0.10608 

5976.53 

Fe II 2344.2140 

0.11400 

0.37211 

5978.53 

Fe II**** 2345.0011 

0.15300 

0.10861 

6016.33 

Fe II*** 2359.8278 

0.06790 

0.07552 

6030.93 

Fe II* 2365.5518 

0.04950 

0.13941 

6053.64 

Fe II 2374.4612 

0.03130 

0.29999 

6074.81 

Fe II 2382.7650 

0.32000 

0.41259 

6091.62 

Fe II* 2389.3582 

0.08250 

0.19033 

6131.92 

Fe II*** 2405.1638 

0.02600 

0.05041 

6133.08 

Fe II** 2405.6186 

0.23700 

0.25156 

6147.44 

Fe II*** 2411.2533 

0.21000 

0.18791 

6148.84 

Fe II**** 2411.8023 

0.21000 

0.10497 

6154.56 

Fe ip*** 2414.0450 

0.17500 

0.10723 

6629.09 

Fe II 2600.1730 

0.22390 

0.44168 

6648.70 

Fe II** 2607.8664 

0.11800 

0.20213 

6660.91 

Fe II* 2612.6542 

0.12600 

0.25575 

6675.56 

Fe II** 2618.3991 

0.05050 

0.10022 

6685.89 

Fe II**** 2622.4518 

0.05600 

0.04794 

6696.08 

Fe II* 2626.4511 

0.04410 

0.15064 

6702.78 

Fe II**** 2629.0777 

0.17300 

0.12269 


Note. —Absorption features of Fe II and Si II ground and fine-structure states 
identified in the spectrum of GRB 051111. The columns are (left to right): (i) 
Observed wavelength, (ii) line identification, (iii) oscillator strength (Prochaska 
et al. 2003), (iv) redshift of the line, and (v) rest-frame equivalent width. 
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Table 2 


Column Densities of Fe II and 


Si II Ground and Fine-Structure States 


Ion 

log IV 

log (N J /N) 

Si II 

> 16.18 


Si II* 

14.96 ± 0.10 

< -1.22 

Fell 

15.21 ±0.10 


Fell* 

13.89 ±0.14 

-1.32 

Fe II** 

13.67 ±0.17 

-1.54 

Fe II*** 

13.52 ±0.15 

-1.69 

Fe Xf**** 

13.18 ±0.20 

-2.03 


Note. —Ionic column densities of Fe II and Si II ground 
and fine-structure state. The columns are (left to right): 
(i) Ion, (ii) logarithm of the column density, and (iii) loga¬ 
rithm of the ratio of fine-structure to ground state column 
density. 




Line profiles of Fe II transitions plotted as a function of velocity relative to the systemic redshift of z = 1.54948 
















Same as Figure 1 but for Fe II*. 
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Same as Figure 1 but for Fe II**. 
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Fig. 4.— Same as Figure 1 but for Fe II***. 
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Same as Figure 1 but for Fe II****. 
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log [T/K] 


Fig. 6.— Contours of temperature and electron density defined by the ratios of the four Fe II fine-structure states to the 
ground state, as well as the ratio of Si II* to Si II. The physical conditions inferred from the Fe II ratios are T ;> 400 K 
and n e ^ 1.6 x 10 3 cm -3 , with the best-fit requiring T ~ 10 3 K and n e ~ 6 x 10 3 cm' 3 . However, these conditions do not 
agree with the observed column density of Si II* from which we derive a much lower density. The thickness of the lines 
designates a 0.15 dex uncertainty in the ratio. 





